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Abstract 



o 

' The Euclidean algorithm makes possible a simple but powerful gen- 

^S) . eralization of Taylor's theorem. Instead of expanding a function in a 

series around a single point, one spreads out the spectrum to include any 
number of points with given multiplicities. Taken together with a simple 
• expression for the remainder, this theorem becomes a powerful tool for 

approximation and interpolation in numerical analysis. We also have a 
corresponding theorem for rational approximation. 



Generalize Rolle's Theorem 



< 

' Let h{x) = ni=i(^~^j)'"' foi' distinct Xi € [a, b] C IR with multiplicity rrii > 1, 

. and let n = deg{h{x)). Given two functions /(x) and g{x), we say that f{x) = 

' h 

d , g(x) mod{h{x)) or f{x) — g{x) if for each 1 < i < r and Q < k < mi 

f''\x,)=g('^\xi). (1) 

If f{x) and g{x) are polynomials, then f{x) = g{x) is equivalent to saying that 
^ ' if /(x) and g{x) are divided by the polynomial h{x) (the Euclidean algorithm), 

CO , they give the same remainder. We denote the factor ring of polynomials modulo 

' h{x) over the real numbers M by 

O ; lR[x]h := IR[x]/ < h{x) >, 

2 ! see [H p.266]. 

I Generalized Rolle's Theorem: Let f{x) G C[a,5] and [n — l)-times 

differentiable on (a, 6). If f{x) — mod{h{x)) , then there exist a c G (a, b) such 
that /(""i)(c) = 0. 

I / 1 u <C 

rN , Proof: Following p. 381, define the function a(u,v) I „' 

?H o ^ t- ,1 ^- ' \ 0, u > V 



The function a is needed to count the simple zeros of the polynomial h{x) and 
its derivatives. 

Let 4f^h'^^^ denote the number of simple zeros that the polynomial equation 
/i('=)(a;) = has. Clearly i^h = r = E[=i cr(0, m,). By the Classical Rolle's 
theorem 

r r r 

Continuing this process, we find that 

r ^ ^ ^ 

4^h" = ^^(2, m,) + {#h') -l=Y ^(2, m,) + ^ m,) + ^ a(0, m,) - 2, 

i—1 i—1 i—1 i—1 

and more generally that 

r r k 

ifh^"^ = Y '^{k, m,) + - 1 = E E ™») - ^ 
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for all integers fc > 0. 

For A: = n — 1 , we have 



^Y.J2 "^0 - (n - 1) = - (n - 1) = 1. 

i=l j=o i=l 

The proof is completed by noting that =fff^''^ > =ffh^''^ for each fc > 0, and hence 

^f(n-l) > 1. 

□ 

Approximation Theorems: 

We can now prove 

Generalized Taylor's Theorem: Let f{x) G C[a, b] and n times differen- 
tiable on (a, b). Suppose that f{x) = g(x) mod{h{x)) for some polynomial g(x) 
where deg(g) < deg{h). Then for every x G [a, 6] there exist a c G (a, 5) such 
that 



f{x)^g{x) + ^^hix). 



n 



Proof: For a given x G [a, b], define the function 

= (2) 

In the case that x — Xi for some 1 < i < r, it is shown below that p{t) has a 
removable singularity and can be redefined accordingly. Noting that 

p{t) = mod{h{t)) and p{x) = 

it follows that p{t) — mod{h{t)(t — x)). Applying the Generalized RoUe's 
Theorem to p{t), there exists a c G {a,b) such that p("^(c) = 0. Using ([2]), we 
calculate ^'■"•'(i), getting 



' h{x) ' 
so that 

= p(")(c) = /(")(c)-[/(:r)-g(:r)]-^ 

/(,(a;) 

from which the result follows. 

□ 

Applying the theorem to the case when x = Xi, we find by repeated application 
of L'Hospital's rule that 

lini /(^)-g(^) ^ /("'H^O - .g^^-Ha^O ^ /^"Hc) 

There remains the question of how do we calculate the polynomial g{x) 

with the property that f{x) = g{x) where deg{g{x)) < deg{h{x))l The brute 
force method is to impose the conditions (H)), and solve the resulting system 
of linear equations for the unique solution known as the osculating polynomial 
approximation to f{x), see [3], [4j p. 52]. A far more powerful method is to 
make use of the special algebraic properties of the spectral basis of the factor 
ring lR[x]h, as has been explained in [S1[S]. See also [7]. 

Much time is devoted to explaining the properties of Lagrange, Hermite, 
and other types of interpolating polynomials in numerical analysis. In teaching 
this subject, the author has discovered that many of the formulas and theorems 
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follow directly from the above theorem. For rational approximation, we have 
the following refinement: 

Rational Approximation Theorem: Let f{x) G C[a, b] and n times 
differentiate on (a, b). Let u{x) and v{x) be polynomials such that v{0) = 1 and 
det(u(x)v{x)) < deg{h{x)), and suppose that f{x)v{x) — u(x) = modihix)). 
Then 

f{x) = ^ + 

v[x) n\v[x) 

for some c G (a, 6). 

Proof: Define the function 

p{t) ^ f{t)v{t) - u{t) - [f{x)v{x) - uix)]^ 

where x G [a,b]. Clearly, p{t) = mod{h{t){t — x)). Applying the Generalized 
Rolle's Theorem to p{t), it follows that there exist a c G (a, b) such that 

f{x)v{x) ~ u{x) = -(-) [f{t)v{t)]t_,M^), 
from which it follows that 

fix) = ^ + -^[f{t)v{t)](-\c)h{x). 

□ 
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